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ABSTRACT 


A  new  methodology  has  been  developed  to  predict  the 
impulsive  noise  generated  by  a  transonic  rotor  blade.  The 
formulation  uses  a  full/potential  finite-difference  method  to 
obtain  the  pressure  field  close  to  the  blade.  A  Kirchhoff 
integral  formulation  is  then  used  to  extend  these  finite- 
difference  results  into  the  farlfield.  This  Kirchhoff  foumula 
is  written  in  a  blade- fixed  coordinate  system.  It  requires 
initial  data  across  a  plane  at  the  sonic  radius.  This  data 
is  provided  by  the  finite/difference  solution.  Acoustic  pres¬ 
sure  predictions  show  excellent  agreement  with  hover  exper¬ 
imental  data  for  two  hover  cases  of  0.88  and  0.00  tip  Mach 
number.  The  latter  of  which  has  delocalised  transonic  flow. 
These  results  represent  the  first  successful  prediction  tech¬ 
nique  for  peak  pressure  amplitudes  using  a  computational 
code.  ^ - 


LIST  OF  SYMBOLS 


=  atmospheric  speed  of  sound,  m/a 
=  rotor  chord  length,  m 
=  nondimensional  time  step 
=  identity  matrix 


=  Jacobian  of  coordinate  transform  matrix 


=  hover  tip  Mach  number 


=  local  Mach  number 


=  absolute  pressure,  N/m2 
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=  nondimensional  scaled  disturbance  pressure, 

(P  -  Po)/(Po«o<V*) 

=  atmospheric  pressure,  N/m2 
=  dimensional  distance  along  the  rotor  in  the 
spanwise  direction,  m 
=  radial  distance  to  the  rotor  tip,  m 
=  free  stream  subtraction  term  in  Eqn.  (5) 

=  nondimensional  time  in  nontransformed  coordi¬ 
nate  system 
=  dimensional  time,  a 
—  scaled  azimuthal  angle  location, 

(l/*)(P  +  0-tan-‘/?) 

=  contravariant  velocity  normal  to  the  17,  f  plane 
=  contravariant  velocity  normal  to  the  £ ,  f  plane 
=  nondimensional  scaled  radial  velocity  compo¬ 
nent,  tr/uiRt 

=  contravariant  velocity  normal  to  the  plane 
=  nondimensional  coordinate  system  aligned  with 
the  blade 

=  blade  radial  location,  -  1]1/3 

=  p'-VJ 

=  inverse  of  the  blade  aspect  ratio,  c/ R 
=  specific  heat  ratio 

=  jump  in  potential  across  the  wake  normalized 
by  (ooc) 

=  blade  angular  speed,  rod/ a 
=  dimensional  velocity  potential,  m2/a 
=  non-dimensional  velocity  potential,  d/cao 
=  fluid  density  normalised  by  free-stream  values 
=  p/J 

=  atmospheric  density,  kg/m * 
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a  -  nondimens  ional  integration  variable,  — x 

r  =  normalized  time  in  the  transformed  coordinate 

system 

8  =  azimuthal  angle  in  cylindrical  coordinates  rela¬ 

tive  to  the  blade 

.  =  transformed  coordinate  system 
INTRODUCTION 

« 

Delocalisation  occurs  for  a  high-speed  rotor  when  the 
supersonic  flow  region  at  the  blade  tip  connects  to  the 
farfield  beyond  the  sonic  cylinder.  The  sonic  cylinder  is  the 
radial  location  beyond  which  the  undisturbed  flow  appears 
to  be  supersonic  to  an  observer  on  the  blade.  Flowfield  non- 
linearities  are  very  significant  at  delocalisation1  and  these 
nonlinearities  cause  large  changes  in  the  far-field  acoustic 
levels. 

Linear  methods  with  monopoles  and  dipoles  have  con¬ 
sistently  underpredicted  the  far-field  acoustic  levels  for  high 
tip  Mach  numbers.  The  methods  developed  by  Ffowcs 
Williams  and  Hawkings2  improve  upon  linear  formulations 
by  including  the  nonlinear  quadrupole  term  in  the  acoustic 
pressure.  Several  methods  have  been  developed  for  evalu¬ 
ating  this  quadrupole  term3-1,  but  difficulties  concerning 
both  the  size  of  the  quadrupole  volume  and  the  acquisition 
of  aerodynamic  data  make  this  term  difficult  to  integrate 
accurately.  To  date,  neither  computer  codes  nor  experi¬ 
mental  techniques  have  been  able  to  successfully  predict 
the  far-field  acoustic  pressures  associated  with  high-speed, 
delocalized  rotor  flows. 

Computational  fluid  dynamics  (CFD)  computer  codes 
in  their  non-appracimated  forms  properly  handle  all  flow- 
field  non-linearities.  Several  computer  codes  exist  for  solv¬ 
ing  the  three-dimensional  full-potential  equation  for  heli¬ 
copter  rotor  flowfields.1-10  Typically  these  codes  are  well 
developed  for  surface  pressure  and  load  predictions  close 
to  the  blade.  Far  from  the  blade  however,  the  grid  den¬ 
sities  are  insufficient  to  resolve  the  details  of  the  acoustic 
pressure  field.  Extending  these  codes  to  predict  acoustics 
is  too  costly  and  perhaps  impossible  due  to  the  huge  grids 
required  and  instability  problems  that  result  from  increas¬ 
ing  the  size  of  the  computational  domain.  These  instabil¬ 
ities  occur  as  the  radial  distance  away  from  the  blade  is 
increased,  thereby  increasing  the  freestream  Mach  number 
relative  to  the  blade.11 


FPR,  which  stands  for  Full-Potential  Rotor  code.  It  solves 
the  fully  conservative,  three-dimensional,  unsteady  form  of 
the  full-potential  eqauation  and  is  described  in  Refs.  8  and 
9.  The  FPR  code  is  used  to  predict  the  pressure  field  both 
on  the  blade  and  along  a  plane  located  at  the  sonic  radius. 
A  computed  solution  at  the  sonic  radius  is  well  within  the 
limits  imposed  by  the  maximum  grid  size  and  instability 
problems. 

In  extending  this  result  to  the  far  field,  we  use  a 
Kirchhoff  integral  formulation  that  has  been  developed  by 
Isom13.  This  Kirchhoff  formulation  uses  integrals  of  the 
finite-difference  results  at  the  sonic  cylinder  to  determine 
the  acoustic  pressures  in  the  far  field.  It  treats  the  linear 
sonic  cylinder  as  a  caustic  where  all  the  important  data 
originates.  This  data  is  then  propagated  along  linear  char¬ 
acteristics  to  the  far  field. 

The  Kirchhoff  formula  uses  non-rotating  surface  inte¬ 
grals  of  pressure  and  velocity.  This  non-rotating  control 
surface  should  be  large  enough  to  enclose  the  rotor  and 
the  region  of  nonlinear  flow.  The  linear  sonic  cylinder  suf¬ 
ficiently  encircles  the  nonlinear  tip  region.  The  flow  field 
on  the  sonic  cylinder  comes  from  a  numerical  solution  us¬ 
ing  a  conservative  full  potential  transonic  code.  This  ap¬ 
proach  has  several  practical  advantages.  A  large  volume 
of  quadrupole  integral  is  replaced  by  a  surface  integral 
over  the  sonic  cylinder  surface.  This  control  surface  is  a 
non-rotating  one,  eliminating  mathematical  singularities  of 
the  Ffowcs  Williams  and  Hawkings  equation.  More  impor¬ 
tantly,  the  dose  proximity  of  the  sonic  cylinder  to  the  rotor 
tip  means  that  the  finite-difference  code  only  needs  to  solve 
the  near-field  region  near  the  tip. 

THE  FULL-POTENTIAL  ROTOR  CODE 

The  Full-Potential  Rotor  solves  the  unsteady  full-po¬ 
tential  equation  in  strong  conservation  form.  This  equation 
is  written  below  in  generalised  coordinates. 


a  a  fV.  a  (,v  9 


with  density  given  by: 


(V  ♦  *,)♦,  -  (W  ♦  *)♦,] 
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The  idea  in  this  paper  la  to  use  a  combined  finite- 
difference  code  and  Kirchhoff  formulation  in  order  to  pre¬ 
dict  the  acoustic  far  field.  The  finite-difference  code  used  is 


U ,  V  And  W  are  contrAvnriAnt  velocities  given  by 


If  “  it  +  Ajffj  +  As#,  +  A§#f 
V  =  ijt  +  +  Aj*,  +  Aflfit  (3) 

^  =  f»  +  As#<  +  Ag4,  +  A|t{ 

where 

Ai  =  (i  +  {,‘+ (2 

Aj  =  nl  +  n\  +  nl 

a»=<:+<:+<:  (4) 

Ag  —  (tla  +  +  (ala 

A*  =*  £«f,  +  f|(t  +  (a$a 

A«  *  fa(a  +  Vjffr  +  Tab 

AU  velocities  Are  normAlised  by  «<»,  distAnces  by  the  Airfoil 
chord  length.  And  time  by  the  combine! ion  (c/sgo).  Density 
is  normalised  by  the  freeetream  value 

Equation  (1)  is  solved  by  using  first-order  backward 
differencing  in  time  end  second-order  centre!  differencing 
in  specs.  The  tempore!  density  derivative  is  locally  lin¬ 
earised  about  the  old  time  levels  in  a  manner  that  preserves 
the  conservative  form.  The  resulting  equation  is  approxi¬ 
mately  factored  into  £,  ty.  and  f  operators.  The  steady-state 
version  of  this  equation  can  be  written  as 


[/  +  hU~t(  ±  Dkt(  -  x 

[/  +  U¥*i,  ±£1*-  j^<{(4A*)"<r](*"+l  -  •")  - 

£  [<«(*!)*  +  l,(AT  +  ww  -  R»] 

(») 

where  i(,  and  S(  represent  central-difference  operators 
in  space.  A  steady  state  ADI  relaxation  algorithm  solves 
this  factored  equation.  The  quantity  JM  represents  a  nu¬ 
merical  truncation  error  term  caused  by  incomplete  metric 
cancellation. 

A  sample  grid  for  a  rotor  calculation  is  also  shown  in 
Fig.  1.  An  O-grid  has  bom  chosen  for  the  bask  finite- 
difference  grid  because  of  its  eflkient  use  of  grid  points. 
Rotor  Sows  are  computed  by  sesigning  an  appropriate  ro¬ 
tational  coordinate  velocity  to  such  grid  point.  These  coor¬ 
dinate  velocities  are  given  by  £«,  q,.  and  {,  Am  a  result,  the 


rotor  and  the  attached  finite-difference  grid  move  through 
■till  air  and  the  velocity  potential,  *,  is  equal  to  tero  in  the 
far  field. 

Boundary  conditions  for  the  calculation  are  shown  in 
Fig.  1.  At  the  outer  radial  boundary  of  the  O-grid,  a  non- 
refiection  condition  is  used  to  prevent  the  accumulation  of 
disturbances  there.  This  boundary  condition  is  described 
in  Ref.  0.  On  the  surface  of  the  blade,  a  transpiration  con¬ 
dition  is  used  to  simulate  angle  of  attack  conditions.  Along 
the  inner  boundary  plane,  normal  to  the  rotor,  the  spanwise 
contra  variant  velocity  V,  is  set  equal  to  Vt- 

For  lifting  cases,  the  shed  vortkity  is  specified  as  a 
jump  in  potential,  I*,  imposed  across  the  coordinate  plane 
£  *  0.  This  coordinate  plane  is  approximately  aligned  with 
the  shear  layer  from  the  trailing  edge.  Density  is  assumed 
to  be  continuous  across  the  cut,  and  an  unsteady  transport 
equation  for  the  potential  jump  is  implemented  in  the  wake. 

The  straamwiee  flux  terms  use  upwind  density  bias¬ 
ing  in  regions  of  supercritkal  flow  to  ensure  stability  of  the 
algorithm.  The  density  ia  biased  in  both  f  and  the  £  di¬ 
rections.  The  density  always  biases  upwind  no  matter  how 
the  grid  lines  are  oriented  to  the  freeetream. 

If  the  finite-difference  grid  is  extended  radially  be¬ 
yond  the  sonk  cylinder,  then  the  freeetream  flow  relative  to 
the  blade  k  supereonk  and  the  full-potential  aquation  be¬ 
comes  hyper boUc.  The  FPR  code  uses  central-differencing 
In  space  with  an  upwind-biaeed  artificial  density  in  regions 
of  supereonk  flow.  This  method  has  been  shewn  in  Ref.  11 
to  give  stability  problems  in  regions  where  the  freeetream 
flow  k  supereonk. 

In  order  to  remedy  these  instabilities  in  the  far  field, 
the  coefficients  D  and  E  have  been  added  in  the  original 
FPR  algorithm  to  act  as  stabilisers.  The  I  k  a  finite- 
difference  operator  chosen  so  that  it  always  adds  upwind 
data  to  the  system.  The  parameters  D  and  E  control  how 
much  damping  k  added  in  the  £  and  c  directions  respec¬ 
tively.  The  effect  of  these  coefficients  k  to  increase  the 
diagonal  dominance  in  the  solution  matrices.  These  terms 
approach  sero  as  the  solution  converges  to  a  steady-state 
result.  Thk  modification  to  the  FPR  code  ensures  a  sta- 
bk  convergence  at  outer  radial  stations  for  cases  where  the 
freeetream  Mach  number  k  well  above  one. 


THE  FAR- FIELD  GEOMETRICAL 
ACOUSTICS  FORMULATION 

Numerical  data  far  the  potential  in  a  region  near  the 
blade  tip  can  be  used  in  a  geometrical  acoustics  approxima¬ 
tion  fat  pressure  in  the  far  field.  This  approximation  will 
not  be  valid  however,  near  the  blade  tip  where  the  fiowfield 
is  non-linear  and  rapidly  varying. 

The  geometrical  acoustics  approximation  can  be  based 
on  the  small-disturbance  potential  equation.  One  form  of 
this  equation  is 


(«) 


The  timo-diffsrentiated  second  order  nonlinearity  Q  is  given 


p(j9,T)-po  = 

_  *>*&M**4/t  ((  V,(r,«)-Vr(o,«)  dads 
(3)(6*/*),r  0jJA  T-c  \T-c\i/z 

[(  P(T,  *)  -  P{o,x)  xgn{T-c)dcdx 

38*0  JJA  T-c  |T-o|‘/» 

_  Poeg(-rH)Ad«V»  ff  P*(T,t)  -  P*(o,t)  dods 
JJA  T-c  T-c 

(8) 

The  location  for  the  pressure  evaluation  is  determined 
by  the  values  of  0  and  T,  where  0  is  a  non-dimensional  ra¬ 
dial  locate tion  and  T  is  determined  by  the  asimuthal  angle 
relative  to  the  blade.  The  use  of  Eqn.  (8)  is  restricted  to 
the  plane  of  the  rotor  blade. 

The  surface  integrals  in  Eqn.  (8)  are  located  at  the 
sonic  cylinder,  (r  *  oo/w).  Choosing  this  location  for  ini¬ 
tial  data  greatly  simplifies  the  result.  The  integrals  at  this 
location  are  evaluated  using  the  computed  results  from  the 
FPR  finite-difference  calculation. 


where  t  is  the  ratio  of  specific  huts  and  so  is  the  sound 
speed  in  undisturbed  air.  Note  that  Eqs.  (6)  and  (7)  aasume 
that  the  local  sound  speed  is  constant.  This  ia  not  the  case 
however  with  the  finite-difference  solution  clou  to  the  blade 
surface. 

The  objective  is  to  use  Eqn.  (8)  to  find  the  disturbance 
pressure  in  the  far  field  in  terms  of  nonlinear  initial  data 
on  some  surface  near  the  blade  tip.  This  is  accomplished 
by  using  a  Kirchhoff  solution  to  Eqn.  (8)  that  converts  it 
to  an  integral  equation  with  Q  as  a  source  term.  Next,  a 
stationary  phase  approximation  is  applied  that  reduces  the 
domain  of  all  integrals  to  that  of  a  data  surface  near  the 
tip.  This  initial  data  surface  is  determined  from  the  finite- 
difference  code  that  solves  the  full-potential  equation  ia  the 
region  does  to  the  blade  tip.  The  advantage  of  the  geomet¬ 
rical  acoustics  is  that  the  finite-difference 

computation  need  not  be  continued  from  the  blade  surface 
ell  the  wsy  into  the  acoustic  field. 

With  these  ideee  in  mind,  Eqn.  (6)  can  be  solved  for 
tbs  acoustic  pressure  at  a  point  in  the  far  field.  Details 
of  this  solution  are  quite  lengthy  and  will  be  described  in 
Ref.  13.  The  resulting  formula  for  acoustic  pressure  is  given 
hy 


INTEGRATION  OF  THE 
FINITE-DIFFERENCE  RESULTS 


The  FPR  solution  produces  pressure  data  at  the  sonic 
cylinder  on  an  O-grid  plane  of  data  perpendicular  to  the 
rotor  axis.  This  solution  is  then  interpolated  onto  a  rect¬ 
angular  mash  for  ease  of  integration.  Since  this  plena  does 
not  Us  on  a  constant  radial  station,  a  small  approximation 
is  introduced  hers  which  greatly  simplifies  the  integration. 
The  integrals  then  use  the  rectangular  mash  data  far  in¬ 
put.  A  bi-cubic  spline  smooths  the  data  between  points  on 
the  mash.  A  domain  of  ±3.5  chords  in  both  the  vertical 
and  horiaoutal  directions  provides  sufficient  input  far  the 
integration. 

All  three  integrals  in  Eqn.  (8)  contain  singularities 
when  T  is  equal  to  c.  In  ail  three  cases,  the  aquations 
are  integrated  ia  the  c  direction  up  to  a  small  distance  ±< 
from  the  singularity.  The  contributions  at  the  singularity 
are  computed  separately.  The  functions  inside  tbs  second 
and  third  integrals  can  ha  shown  to  be  symmetric  sad  odd 
at  the  singularity.  Hence  their  integrated  contributions  are 
sera.  The  function  inside  the  first  integral  is  even,  with  a 
non-canceling  contribution  near  the  singularity.  The  value 
of  this  integral  approaches  the  following  limit  at  the  singu¬ 
larity. 


-4- 


*•/.  r  * 

J-oo 

Thi»  singularity  is  not  insignificant.  It  is  ususlly  the 
sum  order  ss  the  rest  of  the  integral  bat  of  opposite  sign. 
Near  the  pressure  peek  however,  it  suddenly  switches  sign 
end  reinforces  the  rest  of  the  integral  veins.  At  the  pressure 
peek  it  contributes  ebout  25%  to  the  first  integral. 

Awey  from  the  pressure  peek,  the  lest  two  integrals  ere 
ebout  the  seme  megnitude  es  the  first.  At  the  peek  bow* 
ever,  the  first  two  integrals  msy  double  in  their  respective 
values.  The  third  integral,  however,  in  creeses  its  value  by 
almost  two  orders  of  megnitude,  showing  haw  important  it 
is  to  the  overall  solution. 

PREDICTION  OF  IMPULSIVE  ROTOR  NOISE 

Our  rotor-noise  prediction  scheme  it  demonstrated  tor 
two  of  the  experimental  test  cases  presented  in  Ref.  IS.  In 
this  experiment,  e  1/7  scale  model  of  a  UH-lH  mein  rotor 
with  a  N  ACA  0012  airfoil  section  with  a  root-to-tip  washout 
of  10.0*.  This  model  wee  run  et  high  speed  hover  with 
thrust  set  to  e  low  level  in  order  to  minimise  the  woke  in¬ 
fluence.  Acoustic  pressure  dote  wee  measured  in  the  plane 
of  the  rotor  et  a  redial  location  equal  to  SR.  A  variety 
of  hover  tip  Mach  numbers  were  tested  from  M  ■  J  Is 
M  =  1.0 

A  top  view  of  the  finite-difference  grid  for  the  At  ■ 
.0  computation  is  shown  in  Fig.  2.  This  figure  shews  e 
two-dimensional  slice  of  the  grid  in  the  plane  of  the  rotor. 
The  blade  tip  is  located  at  13.7  chords  in  the  spanwiee 
direction.  Only  the  outer  three  chords  of  the  Made  are 
computed  in  the  FPR  solution.  A  large  number  of  grid 
planes  are  clustered  in  the  far  field  in  order  to  obtain  good 
resolution  of  the  acoustic  pressures  at  the  soak  cylinder  (r 
=  15.22  chords).  A  total  of  90  grid  points  were  used  in  the 
chordwiae  direction,  39  in  the  spanwiee  direction,  and  17  in 
the  direction  normal  to  the  surface  of  the  blade. 

An  untwisted,  nonlifting  NACA  0012  rotor  Made  was 
used  for  the  acouetk  predict  too.  This  simplifies  the  finite- 
difference  cakulation  become  the  rotor  woke  can  be  ig¬ 
nored.  Baarweil  at  el.1*  have  shown  that  thrust  level  is 
unimportant  for  high  speed  impulsive  noise  esses.  In  gen¬ 
eral  however,  the  FPR  code  con  incorporate  any  of  a  num¬ 
ber  of  rotor  wake  models. 

Approximately  700  iteration  steps  were  necessary  to 
obtain  a  converged  steady-state  solution  for  this  case.  This 


represents  approximately  10  minutes  of  CPU  time  on  a 
Cray  XMP  4/8  computer.  The  addition  of  the  far-fieid 
damping  terms  causes  this  abnormally  long  run  time.  These 
damping  terms  are  needed  to  control  instabilities  in  regions 
where  the  freestream  is  supersonic.  Typical  iteration  counts 
for  the  FPR  code  on  conventional  grids  are  about  five  times 
lower. 

Computed  Mach  contours  in  the  plane  of  the  rotor  for 
the  Af  =  .9  cose  are  shown  in  Fig.  3.  The  local  sonic 
Mach  number  contour,  Afj  =  1,  deforms  due  to  the  local 
perturbation  velocity  field  and  touches  the  blade  tip.  The 
strong  local  shock  on  the  tip  portion  of  the  rotor  blade  is 
now  extended  beyond  the  blade  tip  to  the  far  field.  This 
phenomenon  is  called  delocalisation.  Delocalisation  causes 
a  stong  shock  in  far-field  acoustic  signatures  at  tip  Mach 
numbers  beyond  the  critical  delocalisation  Mach  number. 
The  critical  delocalisation  Mach  number  for  this  particular 
blade  geometry  is  0.9. 

Figs.  4a  and  4c  shew  computed  premure  distributions 
from  the  FPR  code  on  the  rotor  plane  at  the  sonic  cylinder 
for  the  M  m  .88  and  M  »  .9  cases  respectively.  Only  the 
in-plane  pressure  distributions  are  shown  in  these  plots  al¬ 
though  the  solution  ie  actually  obtained  sc  roes  a  plane  at 
the  sank  cylinder.  These  finite-difference  predktions  are 
then  need  ee  initial  data  far  the  Kirchoff  integral  formula¬ 
tion  in  Eq.  (8)  and  the  solution  k  carried  to  the  far  field. 

The  integrations  in  Eq.  (8)  require  initial  data  over  the 
entire  souk  cylinder.  The  FPR  solutions  only  provide  this 
data  ±3.5  chords  along  the  sonic  cylinder  in  both  the  x  and 
the  ■  direction.  Outside  af  thk  region,  the  pressure*  and 
radial  velocities  are  assumed  to  be  aero.  Thk  k  clearly  the 
ceee  in  the  chordwiae  direction  os  shown  in  Fig*.  4a  and  4c. 
Her*  the  pressure  k  seen  to  fall  quickly  to  sero  inside  the 
chardwks  computatinal  domain  for  the  FPR  cod*.  Peak 
pressure  along  the  vertical  axk  along  the  sonic  cylinder  also 
folk  rapidly  to  freest  r*sm  values.  Thk  k  demonstrated  in 
Fig.  5  where  the  relative  pressure  has  dropped  to  sero  well 
within  the  ±3.5  chord  domain  af  the  FPR  cod*. 

Figs.  4b  and  4d  compare  computed  far-field  pressure 
results  to  experimental  microphone  measurements  at  the 
r  w  SR  radial  station.  Good  agreement  k  seen  between  the 
predktion  and  the  dots  for  peak  pressure  magnitude*  and 
pressure  pule*  waveforms.  For  the  M  *  .88  case  (Fig.  4b) 
the  computed  result  underpredicts  the  peak  amplitude  and 
shows  a  slightly  wider  pule*  width  compared  to  the  exper¬ 
imental  data.  The  symmetrical  pulse  shape  agrees  quite 
well  with  the  data  however. 
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The  experimental  data  for  the  Ai  *■  .9  caee  (Fig.  4d) 
show  a  strong  shock  in  the  pressure  data  that  is  not  seen 
in  the  prediction.  The  prediction  does  show  a  steeper  slope 
on  the  compression  side  however.  This  discrepancy  in  the 
shock  prediction  may  be  attributed  to  two  causes.  First,  the 
Kircboff  integral  formulation  has  an  assumption  of  constant 
sound  propagation  speed  as  shown  in  Eq.  (8).  The  far- 
field  nonlinear  phenomenon  may  steepen  the  waveform  as 
described  by  bom14.  If  this  is  the  case,  then  the  Kirchoff 
integral  formulation  cannot  be  expected  to  model  this. 

A  second  cause  for  disagreement  in  the  shock  structure 
may  be  that  the  FPR  finite-difference  solution  has  smeared 
an  existing  shock  at  the  sonic  cylinder.  Inadequate  grid 
resolution  could  be  the  cause  for  this.  If  this  is  the  ease, 
then  the  resulting  far-field  pressure  pulse  will  also  show  a 
smeared  shock.  In  spite  of  these  problems,  the  good  agree¬ 
ment  between  experiment  and  computation  appears  to  val¬ 
idate  our  new  model  for  far-field  rotor  acoustic  prediction 
techniques. 


CONCLUSIONS 

A  new  method  has  been  developed  to  predict  the  im¬ 
pulsive  noise  signature  generated  by  a  high-epeed  transonic 
rotor  blade.  The  matching  of  a  finite-difference  method 
and  a  Kirchhoff  integral  formula  combines  the  advantages 
of  each  formulation  in  the  regions  where  eech  is  most  ap¬ 
propriate.  The  finite-difference  solution  is  used  close  to  the 
blade,  where  Bow  held  nonlinearities  are  large.  The  Kirch¬ 
hoff  integral  is  used  in  the  far  field  where  its  simplicity 
allows  its  application  over  a  large  domain. 

Predictions  with  this  method  have  shown  excellent 
agreement  with  experimental  data  for  two  hover  casee.  One 
of  these  casee  exhibits  strongly  delocalised  Bow  These  re¬ 
sults  represent  the  first  successful  prediction  technique  for 
peak  pressure  amplitudes  with  a  computational  code. 
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Z  DISTANCE  (CHORDS) 

Fig.  5  Preaaure  prediction  in  the  s  direction.  Prewar*  location  is  on  the  sonic  cylinder  at  midchord. 
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